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Stress transmission in planar open-cell cellular solids is 
analysed using a recent theory developed for marginally rigid 
granular assemblies. This is made possible by constructing a 
one-to-one mapping between the two systems. General triva- 
lent networks are mapped onto assemblies of rough grains, 
while networks where Plateau rules are observed, are mapped 
onto assemblies of smooth grains. The constitutive part of 
the stress transmission equations couples the stress directly 
to the local rotational disorder of the cellular structure via a 
new fabric tensor. An intriguing consequence of the analysis 
is that the stress field can be determined in terms of the mi- 
crostructure alone independent of stress-strain information. 
This redefines the problem of structure-property relationship 
in these materials and poses questions on the relations be- 
tween this formalism and elasticity theory. The deviation of 
the stress transmission equations from those of conventional 
solids has been interpreted in the context of granular assem- 
blies as a new state of solid matter and the relevance of this 
interpretation to the state of matter of cellular solids is dis- 
cussed. 

62.20.-x, 61.43.Gt, 81.40.Jj 
I. Introduction 

Planar foams consist of partitions of space by irregular 
cells and are observed in many natural and man made 
systems. Their ubiquity in nature is often due to evo- 
lutional processes that have led to an optimization of 
mechanical or transport properties against weight con- 
straints. Much effort has been invested in attempts to 
follow the wisdom of Mother Nature and develop materi- 
als with biologically inspired structures. This effort has 
resulted in a growing body of knowledge on the charac- 
terization of the structure as well as on understanding 
of the geometry of the cells in terms of the nonequilib- 
rium dynamical processes that generate them. Essential 
to this entire field is a basic understanding of the rela- 
tions between the internal structural characteristics of 
such media and their macroscopic properties [1]. Such 
an understanding would be of use to a wide range of ap- 
plications, such as design of high-perfrmance and smart 
materials, bio-artificial technology, and the food insus- 
try, to name only a few. A corner stone in the study 
of cellular materials is the work of Gibson and Ashby [2] 
who employed various approximations on ordered cellular 
structures to estimate the mechanical properties of cellu- 
lar networks in general. While that work has significantly 
improved the theoretical basis for stress calculations, it 
has failed to provide a rigorous relation between the lo- 
cal microstructure in general disordered systems and the 
local stress field. It has become increasingly evident in 



recent years that the lack of such a first principle relation 
stands in the way of progress both on the theoretical and 
the technological fronts. 

It is conventional to relate constitutive structural data 
in solids to macroscopic properties by integration over 
microscopic variables [3] [4] [5] [6] [7]. In this approach 
one starts from a given constitutive stress-strain rela- 
tion on a judiciously chosen microscopic lengthscalc and 
coarse-grains this relation until it applies to a suitable 
macroscopic lengthscale. This paradigmatic approach 
has been automatically adopted in recent era for cel- 
lular solids but the main problem has always been the 
translation of the local cellular structure, which can be 
arbitrary in general systems, into local elastic properties. 
The question whether it is really necessary to go down 
this route at all has not been seriously addressed. The 
results in this paper suggest an alternative approach that 
obviates the problem of finding the elastic constants al- 
together. 

On the conceptual level, the problem is the following: 
Suppose a cellular system, whose structure is known in 
every detail, is in mechanical equilibrium under a set of 
forces g acting on its boundaries. Suppose also that the 
properties of the material constituting the network are 
well characterized. The questions are: (i) how to describe 
the discrete field of forces transmitted by the cell walls 
as a continuous stress field? (ii) what are the equations 
that govern this field (the stress transmission equations)? 
and (iii) how do different microstructures affect the local 
stress? In d dimensions it takes d 2 conditions to deter- 
mine the d x d stress tensor a. Force balance gives d 
conditions, 

W = g (1) 

and blanace of torque provides further d(d— l)/2 condi- 
tions - one for each axis of rotation, 

°ij = O-ji (2) 

Equations (1) and (2) give altogether d(d + l)/2 rela- 
tions and therefore additional d(d — l)/2 equations (one 
in d = 2 and three in d = 3) are required to solve for the 
stress components. These 'missing' equations must con- 
sist of information about the constitution of the material 
- the constitutive equations. In conventional elasticity 
theory, this second set of equations is obtained by im- 
posing compatibility conditions on the strain field and 
relating the stress to the strain through constitutive re- 
lations that are independent of the balance conditions. 
While equations (1) and (2) couple the stress to the ex- 
ternal loading, the constitutive equations relate it to the 



properties of the medium and determine how the stress 
fluctuates due to spatial variations of these properties. In 
cellular systems, these variations are governed by the mi- 
crostructure and so it seems that there is no escape from 
the problem of structure-property relations. Yet, what is 
the property that the structure needs relating to? The 
basic issue is not how to relate the structure to clastic 
constants, which we can then plug into the comaptibility 
relations. Rather, it would suffice to find d(d — l)/2 re- 
lations that couple the structural information directly to 
the stress field. 

A development in the study of granular systems has 
provided an unexpected key for progress on this problem. 
In a recent paper, Ball and Blumenfcld (BB) [8] have 
formulated a new theory for stress transmission in two 
dimensional marginally rigid granular assemblies [9]. In 
particular, they have derived a constitutive equation that 
couples the local stress to the local geometry. Marginal 
rigidity is a new state of solid matter, whose existence 
has been verified experimentally [10]. In this state the 
number of intergranular contacts in the system allows to 
determine the microscopic force field from statics alone 
and therefore the macroscopic stress is independent of the 
stress-strain data. BB have further shown that the new 
theory predicts that the symmetric divergence-free stress 
field is the Airy stress function [11], recovering a result 
of conventional elasticity theory and so confirming the 
redundancy of the stress-strain relations in marginally 
rigid systems. 

The aim of this paper is to to adapt the BB formalism 
to planar trivalcnt cellular structures and derive a con- 
stitutive relation between the microstructure of a model 
of such systems and the stress field that develops in them 
under external loading. To this end, an exact mapping 
is formulated between cellular networks and granular as- 
semblies at the marginal rigidity state. By decoupling the 
stress from the strain, the new theory obviates the prob- 
lem of structure-property relationship in the traditional 
sense and provides a direct structure-field relationship. 
In the context of granular systems, marginal rigidity is 
a new state of solid matter and this begs the question 
whether cellular solids and solid foams are also in this 
state. 

The paper is structured as follows: Section II presents 
an argument that leads to the conclusion that the stress 
field in cellular networks at equilibrium can be deter- 
mined from the loading and the microstructure alone, 
independent of any strain-stress information. Section III 
formulates a one-to-one mapping between planar triva- 
lent cellular networks and marginally rigid assemblies of 
grains. Section IV describes the derivation of the con- 
stitutive equation and section V concludes with a brief 
discussion of the results and some of their implications. 

II. The stress field can be determined indepen- 
dently of stress-strain relations 

In this section it is shown from first principles that the 




FIG. 1. The mapping from a trigonal cellular network onto 
an assembly of grains under external loading. The original 
cellular network is shown in dashed lines. The midpoints of 
the cell walls are joined to form triangles around each vertex 
(black circle) . The grains of the equivalent granular assembly 
are centered on the vertices with the corners of the triangles 
representing contacts between neighboring grains. Two such 
virtual grains in contact are shown. 

stress field can be determined from statics alone. The 
argument is slightly different for general trivalent net- 
works, where the cell walls meet at arbitrary angles, and 
for Plateau networks. Therefore these two cases will be 
treated separately throughout the paper. In the follow- 
ing we model a cellular solid by a collection of vertices 
at given positions in the plane, which are connected by 
cell walls whose thickness is much smaller than the size 
of a cell (a dry foam). These cell walls may be straight 
lines or curved as a consequence of the process that has 
given rise to the microstructure and with which we are 
not concerned in the present discussion. 

(a) General networks 

Consider a planar open-cell trivalent cellular network 
that is at mechanical equilibrium under a set of external 
forces, g, acting on its boundaries and which consists of N 
cells. The loading gives rise to tension (or compression, 
which will be regarded in the following simply as negative 
tension) in the cell walls and to torques on the vertices 
of the network. The former are balanced by forces along 
the cell walls, while the latter must be balanced by forces 
that are normal to the walls. By Euler's theorem [12] for 
iV > 1 a cell has on average six neighbors and, since 
each cell wall is shared by two cells, there are altogether 
3N cell walls. Each vertex in the network is the meeting 
point of three walls and each wall connects two vertices 
and so there are exactly 2N vertices. This enumeration 
neglects boundary effects, which introduce corrections of 
order 0(1/N), but which can be readily included into 
the argument that follows and therefore do not affect its 
generality. 

Let us draw around each vertex a triangle whose cor- 
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ners are located at the midpoints of the walls that branch 
from it (e.g. triangle ABC in figure 1) and let us imagine 
that the edges of these triangles are made of rigid struts. 
Each of the 2N triangles is in contact with three neigh- 
bors and the contact points are considered rigid so that 
the triangles cannot rotate relative to each other around 
them. Two triangles in contact exert a force on each 
other whose two components are defined to be the ten- 
sion along the wall and the force normal to the wall that 
contributes to balancing the torques on the vertices due 
to the external loading. Within this description there 
are no torques on the contact points (which may result, 
for example, from a finite thickness of the cell walls). 
Thus defined, the triangles form a rigid frame that trans- 
mits the same forces and moments as the original cellular 
structure. There are 3N contact forces and therefore, to 
determine the intercellular forces, one needs to solve for 
6N unknowns altogether. 

How many equations are available to determine these? 
The triangles are at mechanical equilibrium and so each 
must obey two equations of force balance and one equa- 
tion of torque balance. Multiplied by the 2N vertices, 
this gives 6N equations, exactly enough to determine all 
the forces. The balance equations depend only on the de- 
tails of the microstructure and on the external loading. 
It follows that the intercellular forces can be uniquely 
determined from statics alone. 

(b) Plateau networks 

In Plateau trivalent networks the cell walls meet at 
each vertex at angles of 2tt/3, a structure that abounds 
in nature and man-made materials. In contrast to the 
general case, force balance at every vertex dictates here 
that all the cell walls are under the same tension. This 
leaves the normal force on each cell wall as the only un- 
known and the number of degrees of freedom is reduced 
from 6N to 3N. We use the same construction as in fig- 
ure 1, but instead of considering the frame of triangles, 
let us regard it as consisting of a set of connected rigid 
polygons, each enclosed inside a cell. These polygons 
connect to each other at the same points that the trian- 
gles do, namely, at the midpoints of the cell walls. If we 
define the forces transmitted between every two neigh- 
boring polygons as the forces normal to the cell walls at 
the midpoints, then this rigid frame supports the same 
stresses as the original system. Each of the TV polygons 
has to obey three balance equations, two of force and one 
of torque, which gives altogether 3N equations, exactly 
enouch to determine all the forces. 

So, the intercellular forces in both general and Plateau 
networks can be derived from statics alone and there- 
fore these are isostatic systems. Now, the macroscopic 
stress is just a coarse-grained continuous description of 
the discrete field of forces on the cellular scale. It fol- 
lows that since the latter can be determined without ref- 
erence to compliance then so can the continuous stress 
field. It is evident then that the second set of equations 
that complements the balance conditions must arise from 



the structural characteristics alone. This observation 
has a significant ramification: regardless of the precise 
form that the consitutive equations take, the fact that 
the stress field equations involve no information on the 
stress-strain response means that the stress transmission 
in cellular solids is governed by physical laws that are dif- 
ferent than in conventional elastic solids [13]. This basic 
realization appears to have been neglected in the large 
body of literature on solid foams. Most of the work in 
the field employs stress-strain relations, thus linking the 
stress to the displacement that occurs as the system ap- 
proaches mechanical equilibrium. The fundamental issue 
underlying the present analysis is that in planar cellular 
solids the stress can be decoupled from any deformation. 
Thus, after all displacement has stopped, the stress can 
be calculated from the microstructure and the boundary 
conditions alone. 

III. Mapping cellular networks onto granular as- 
semblies 

A planar trivalent foam can be characterized by the 
positions of all its vertices, i>i, a connectivity matrix that 
specifies which vertices are connected and the curvatures 
of the connecting walls. The cell walls tesselate the sys- 
tem into a set of cells, Cj. It is useful again to discuss sep- 
arately the mappings of general and Plateau networks. 

(a) General networks 

Let us consider again the construction shown in fig- 
ure 1 and note that, although the cell walls may be 
curved, the edges of the constructed triangles are always 
straight. The above enumeration of equations versus de- 
grees of freedom implies that one needs 2N rigid bodies 
and therefore suggests to center the equivalent grains on 
the vertices. The contacts between neighbor grains are at 
the midpoints on the cell walls between them where the 
respective triangles touch. The force f vv > that is trans- 
mitted at the contact between the grains on vertices v and 
v' is defined to be identical to the cell wall at this point. 
Namely, it consists of the tension along the wall, T vv >, 
and the orthogonal component F vv > that assists in bal- 
ancing the torques on vertices v and v 1 [14]. The grains 
are considered to be static under the local forces applied 
on them (i.e., of sufficient rigidity to support the forces 
without deformation), and of sufficient roughness to pre- 
vent slip at any of the contacts under the intergranular 
forces. Note that the sides of triangle v connect the con- 
tact points of the virtual grain centered on vertex v. 

(b) Plateau networks 

Using again the construction of figure 1, let us now 
place a virtual grain inside every cell rather than around 
the vertices. The reason is that this system should pro- 
vide 3N balance equations, which can only come from N 
rigid bodies and N is exactly the number of cells. The 
grains are again chosen to be sufficienly rigid to support 
the forces and the contacts between neighboring grains 
are at the midpoints of the cell walls. At the contact 
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FIG. 2. The vectors f cv and p cv shared between vertex v 
and cell c. f cv connects two neighboring wall midpoints and 
is one edge in a clockwise-directed triangle around vertex v. 
The vector p cv points from the centroid of triangle v to the 
centroid of cell c. 

points the surfaces of the grains are presumed to be nor- 
mal to the 37V unknown forces, and tangent to the wall 
directions. In contrast to the general case, the grains are 
supposed to be perfectly smooth and so incapable of sup- 
porting tangential forces at the contacts. Thus defined, 
the cell wall forces are mapped onto the normal forces 
that neighboring smooth grains exert on each other. The 
mean coordination number is six per grain, each of which 
is then associated with one unknown - the force normal 
to the surface. 

Before we continue it is worthwhile to comment that 
there is a marked difference between real cohesionlcss 
granular assemblies and those described above. Real 
grains can only support compressive forces and tensile 
forces take them out of mechanical equilibrium. In con- 
trast, cell walls can sustain either direction offerees. Cor- 
respondingly, the virtual assemblies described here which 
cellular networks map onto, can also have either type of 
forces between the virtual grains [15]. Although the orig- 
inal BB formalism has been developed for real systems 
under compressive forces, it can readily accommodate an 
overall change of sign into tensile forces. 

IV. The missing constitutive equation 

This section derives the constitutive equation by 
adapting the BB theory to the context of cellular solids. 
The derivation follows three main steps: First, the inter- 
cellular forces are described in terms of a new reduced 
set of forces that satisfies the force balance conditions. 
Second, the new forces are represented in terms of a con- 
tinuous force field. The third and final step is to describe 
the stress as a function of the continuous field and apply 
the requirement of torque balance to obtain the missing 
constitutive equation. 

(a) General networks 



Starting from the construction shown in figure 1, let us 
introduce several definitions. First, we assign directions 
to the sides of the triangles around the vertices so that 
every triangle circulates its vertex anticlockwise. Each 
side is now a vector, r, that lies between a particular 
vertex v and a particular cell c and can be uniquely in- 
dexed r cv (see figure 2). We further define the centroids 
of the cells as the mean positions of the wall midpoints. 
Finally, let us define vectors R cv that extend from the 
centroid of triangle v to the centroid of a neighbor cell 
c. The vectors r cv and R cv form two self-dual networks, 
each edge of which intersects an edge of the other. 

Step 1: Each cell is assigned a force f c that acts at 
the centroid of the cell. It is in terms of the cell forces 
that we wish to parameterize the contact forces. Every 
contact sits between two cells, e.g. the contact point 
between vertices v and v' that is party to cells I and m 
in figure 3. The force that triangle v exerts on triangle 
v' is parameterized as 

fvv' — fl fm (3) 

The sign convention in cq. (3) is such that f c is pos- 
itive when the directed loop formed by the vectors f cv 
around cell c is toward vertex v' at the W-contact and 
vice versa for vertex v. This parameterization achieves 
several aims: (i) By construction, f v / v = f m — fi = —f vv > 
and therefore Newton's third law is satisfied at every con- 
tact, (ii) The forces on each triangle balance automati- 
cally. This can be verified by summing all the forces on 
triangle v', J2v fw- Each of the cell forces /;, f m , and 
/„ contributes to this sum exactly twice, once with a 
positive and once with a negative sign and therefore this 
sum vanishes identically (iii) The number of degrees of 
freedom has been reduced from 6iV to 27V, namely, to 
one force per cell. This coarse-graining corresponds to 
a reduction in the number of equations since the force 
balance constraints have been now satisfied. Only one 
torque balance condition per triangle remains, giving 2iV 
equation to determine the unknowns. 

Step 2: We now wish to continue the discrete cell forces 
into a continuous field, Regarding every cell force as 
located at the centroid of its cell, let us interpolate the 
three forces around every triangle piecewise linearly. The 
interpolation consists of a planar triangular surface whose 
corners are located above the centroids of cells I, m, and 
n, and whose heights are, respectively, f m , and /„. 
The value of the continuous force field at the vertex v, 
is the weighted mean of the cell forces around it, as 
illustrated in figure 4 for the x-component. In terms of 

the cell forces around vertex v can be presented as 

f c = $ v + R C v ■ V$ (4) 

Step 3: In terms of the vectors f cv and forces f c , the 
force moment on triangle v can be written as 

S^J2 r ^c 3 (5) 
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where the sum runs over all the c-cells surrounding grain 
v. Substituting (4) into (5) and observing that ^2 c f cv = 
gives 
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FIG. 3. The parameterisation of the wall forces around ver- 
tex v' in terms of its neighboring cell forces. f vv i = fi — f m is 
the force that triangle v exerts on triangle v'. The sign con- 
vention is that fi is positive because the directed loop formed 
by the vectors r around cell I is toward triangle v' and vice 
versa for cell m. 




FIG. 4. The discrete field of cell forces is made into a con- 
tinuous force field <3> by constructing a piecewise linear tri- 
angular surface around each vertex. An example of such a 
surface for the z-component of the cell forces around vertex 
v is shown. The corners of the triangular surface are above 
the centroids of the cells I, m and n that are located at r ; , 
r m , and f n . The heights of the planar surface at these points 
are, respectively, ff, f^, and f„. The value of the surface 
above the vertex gives In terms of the continuous field, 
the value of any cell force f c , immediately neighboring vertex 
v, is exactly f c = $ + R cv ■ V$. 



Sl j - E r cv R t d ^ j = (<% • V$) 13 (6) 

c,k 

where 

C v — ^ ^ r cv R cv (7) 

c 

and the component notation has been suppressed. C v 
is a fabric tensor that characterizes the local geometry 
around vertex v and plays a crucial role in the following 
analysis. It is convenient to separate C v into its anti- 
symmetric part A v e — (C v — C^)/2 and symmetric part 
P v = (C v +C^)/2, where A v is a prefactor and e = (^/n) 
is the 7r/2 rotation tensor. Both these parts have direct 
geometrical interpretations: A v — \r cv x R cv \/2 is the 
area associated with vertex v that is shown shaded in 
figure 2. The entire system can be conveniently tesse- 
lated and covered by these areas, namely, A sys = J2 V A v - 
This feature is useful for calculations of effective volumes 
[16]. To identify the geometrical interpretation of the 
symmetric part, let us rewrite the vectors f cv and R cv in 
expression (7) in terms of the vectors s cv and t cv shown 
in figure 2. It is straightforward to show that P v can be 
represented in the form 




c 



The departure of P v from zero around vertex v is a direct 
measure of how much the v triangle is rotated relative to 
its immediate environment. For example, if for a partic- 
ular self-dual pair the vectors f and R are perpendicular 
then the symmetric part is exactly zero. Thus, this ten- 
sor characterizes the rotational disorder around vertex v. 
Expression (8) makes it straightforward to see that sum- 
ming P v over any patch of the system the contribution 
from cells fully enclosed inside the patch cancel out and 
the only non-vanishing contribution comes from bound- 
ary terms. Thus, the average of the symmetric part over 
the entire system vanishes at least as fast as 1/ \/N. 

Consider now a small patch T with boundary dT that 
consists of several cells. The mean stress inside this re- 
gion is 

ar = (S v ) r /A T = (C v ■ V$) r /A T (9) 

where (. . .)r stands for an area average over the patch T 
and Ar is its area. It has been shown in [8] that 

(C v ■ V$) r - (C v ) r ■ (V$) 9r (10) 

where (. . .}gr stands for an average over the boundary of 
T. Combining (9) and (10) we have 
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FIG. 5. Making all the grains rough and inserting a vanish- 
ingly small ball bearing (white circle) between touching grains 
at the contact maintains the original intergranular forces. 

(V$) ar - (C v )^A r a r = (11) 

where the last equality uses the fact that the average of 
(Cv)r = Ar[e + (P)r] and that the average (P) vanishes 
faster than 1/y/N. We can now apply the yet unused 
condition of torque balance to a small region in the sys- 
tem: 

A{d v -f$} = A{c v r 1 a} = o (12) 

where A{. . .} stands for the antisymmetric part of {. . .}. 
Separating C v into its symmetric and antisymmetric 
parts we obtain two equations: The first, 

A{&} = 

reassuringly yields the expected condition that the mean 
stress field is symmetric. The second, 

A {Pe^a} = p xx a yy + p vv <j xx - 2p x y<j xy = (13) 

is new and gives the missing constitutive equation. By 
defining a rotated version of P, Q = ePe^ 1 , this equation 
can be rewritten more compactly 

Q:a = (14) 

This new constitutive equation couples the stress to local 
fluctuations in the geometry of the cellular structure and 
completes the set of equations for the stress field in planar 
systems. Just like the tensor P, Q characterizes the local 
rotational disorder and it is this particular aspect of the 
local geometry that the stress field couples to. 
Application to Plateau networks 

As formulated above, the theory applies only to gen- 
eral networks whose cell walls meet at arbitrary angles. 
This is because in general networks there are 6N degrees 



of freedom to determine on the cellular level. This ne- 
cessitates 2N rigid bodies for obtaining equations for the 
balance of forces and torques, which are conveniently pro- 
vided by the triangles around the vertices. In contrast, 
loading a Plateau network gives rise only to 3N degrees 
of freedom on the cellular level: the mean normal force 
on every cell wall. This then implies that we have to use a 
mapping to assemblies of smooth grains that are located 
within the cells and which transmit at each contact only 
force that is normal to the contact surface. The adapta- 
tion to Plateau networks goes through the observation [8] 
that it is possible to extend the theory to assemblies of 
smooth grains. This is done via a further mapping from 
assmeblies of N smooth rigid grains into assemblies of AN 
rough rigid grains as follows: First, the original smooth 
grains are assigned an infinite friction coefficient, making 
them rough and thus, in principle, able to support tan- 
gential contact forces. Second, an infinitesimally small 
circular grain is inserted at each contact point, as shown 
in figure 5. The function of these additional grains is 
to act as ball bearings and so transmit only the normal 
forces between the original grains. This preserves the 
directions and magnitudes of the forces in the original 
assmebly (see figure 5). The mean coordination number 
in the original system is six per grain giving 3N contacts 
altogether. This means that there are 3N extra grains 
and the new assembly consists of AN grains all in all. The 
new system reproduces the original field of intercellular 
forces in the smooth system in the limit when the radii 
of the ball bearings tend to zero. Let us make sure that 
the number of equations still matches the number of un- 
knowns: In the original system there were 3A contacts, 
each of which corresponding to one unknown. In the 
'decorated' system each contact point split into two and 
so we have 6N contacts. Since each ball bearing has two 
contacts then the force balance equations on them reduce 
to the condition that the force at one contact point is the 
same magnitude as the force at the opposite contact, but 
opposite in direction. This immediately satisfies torque 
balance on the ball bearing. There remain therefore the 
3A equations of force and torque balance on the original 
grains, which exactly take care of the 3N unknowns. 

Now we have an assembly of rough grains for which 
the theory is valid and to which equation (13) applies. 
Choosing the radii of all ball bearings to be the same, 
BB have shown that the equation remains non-trivial 
when this radius tends to zero. This extends the theory 
to assemblies of smooth grains and therefore to Plateau 
networks. 

V. Conclusion, discussion and future directions 

This paper has addressed the problem of determining 
the stress field in planar trivalent cellular solids that are 
in mechanical equilibrium under an external loading. By 
comparing the number of degrees of freedom to the num- 
ber of available equations on the cellular level it has been 
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shown that the intercellular forces can be determined 
from statics alone, namely, in terms of only the exter- 
nal loading and knowledge of the microstructure of the 
system. It has been argued then that since the continu- 
ous stress field is only a coarse-grained description of the 
discrete field of intercellular forces and since those can be 
determined independently of material compliance, then 
the stress field must be independent of any stress-strain 
relations. This led to the question what replaces the 
conventional constitutive equation of elasticity theory in 
the determination of the stress field. To find the missing 
equation, a mapping has been constructed between pla- 
nar trivalent cellular networks and granular assemblies 
of infinitely rigid grains at the so-called marginal rigidity 
state. General structures, where cell walls meet at arbi- 
trary angles, map onto assemblies of rough grains that 
are centered on the vertices of the original cellular net- 
work, while Plateau structures, where cell walls meet at 
angles of 27r/3 everywhere, map onto assemblies of per- 
fectly smooth grains that are enclosed inside the original 
cells. The mapping enables the application of the Ball- 
Blumcnfcld theory for stress transmission in assemblies 
of rigid grains to solid foams. In particular, a consti- 
tutive equation has been derived using this approach, 
which couples the stress to the local rotational disor- 
der of the cellular microstructure. This relation replaces 
the conventional stress-strain relations and redefines the 
problem of structure-property relationship in these sys- 
tems. The analysis presented in this paper has several 
fundamental and practical implications, some of which 
are listed below. 
Fundamental issues: 

(i) Although the discussion here focused on comparison 
with elasticity, it should be remembered that all conven- 
tional theories of stress transmission are based on consti- 
tutive relations that link the stress to deformation. Elas- 
ticity employs strain while elastoplastic and viscoelas- 
tic theories are based on relations to strain rate. By 
resorting to statics alone, the theory developed here is 
therefore outside the existing paradigm altogether. It 
has been shown by BB that the new formalism recov- 
ers the Airy stress function as the solution for the sym- 
metric divergence-free stress field. The fact that this 
solution, which is usually obtained within conventional 
elasticity theory, can be derived without the use of any 
stress-strain relation clearly demonstrates that these are 
redundant for determining the stresses that develop in 
trivalent solid foams. This also means that there is a 
connection between the two approaches. To see why such 
a relation must exist suppose we apply an infinitesimally 
small strain to a solid foam, of the type discussed here, at 
the end of which the medium settles into the aforemen- 
tioned isostatic state. Elasticity theory should be able to 
calculate the stress field just before this state is reached, 
while the new theory should be able to do so after all 
deformation has ceased. The stress field cannot change 
discontinuously at the end of the infinitesimal straining 



process and therefore both theories must yield the same 
answer in the limit when the strain goes to zero. This 
issue is currently under investigation by this author and 
preliminary results suggest that indeed there is a relation 
between the two theories. 

(ii) As mentioned already, in the context of granular as- 
semblies, when all the intcrgranular forces can be deter- 
mined from statics alone the system is in a marginal rigid- 
ity state [10] [17], which is interpreted as a new state of 
solid matter. In granular assemblies, the basic condition 
for being in this state is that the average coordination 
number per grain is three for rough grains and six for 
smooth. The state of matter of any system can be de- 
termined by its response to external mechanical loading: 
Gases need to be held together by an external pressure, 
liquids flow under stress and solids first deform and then 
transmit the stress statically. The difference between the 
equations of stress transmission in marginally rigid as- 
semblies and in conventional solids implies different re- 
sponses to external loading and hence, by definition, that 
these are in different states of solid matter. In gran- 
ular systems this new state lies between the fluid and 
the traditional solid states. Adding more contacts be- 
tween grains necessitates additional constitutive informa- 
tion beyond the knowledge of the structure alone, while 
a reduction in the number of contacts takes the assembly 
outside mechanical equilibrium and renders it fluid. This 
feature of existence between two conventional states ap- 
pears to have no immediate analog in solid foams. It is 
possible that such an analog can be found in the context 
of the rheology of foams close to the so-called jamming 
transition point. It is interesting to note that in granular 
systems there is a direct relation between the constitutive 
equation (14) and the dynamics of yield [18]. It would be 
interesting to try to extend this connection to the flow of 
foams. 

(iii) The basic principle of force determination on the 
cellular level without resorting to the compliance of the 
cell walls does not appear to rely heavily on the two- 
dimensional nature of the system. In fact, it is possible 
to extend this argument to open-cell foams where each 
vertex is the meeting point of four edges and six walls, 
a rather ubiquitous class of structures in nature. This 
suggests that the stress field in such three-dimensional 
systems should also be independent of stress-strain rela- 
tions and therefore that there is a larger class of solid 
foams that are at a state of marginal rigidity. 
Practical issues: 

(i) First and foremost, this theory will clearly make it 
possible to model disordered planar cellular materials to 
predict stress transmission, extending the current capa- 
bilities. In particular, the theory should make it straight- 
forward to predict 'hot spots' of stress concentration in 
such systems and make useful connections between the 
occurrence of hot spots and the microstructure. This 
knowledge is expected to both improve prediction of fail- 
ure in advanced materials and set new goals for the de- 
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sign of microstructural characteristics that prevent such 
vulnerability. 

(ii) The formal mapping between granular assemblies and 
cellular systems can be used to transfer results between 
the two fields. Granular assemblies can only approach 
the marginal rigidity state under careful consolidation 
processes [10]. In contrast, it is quite common for the 
dynamics of formation of foams to give rise to triple junc- 
tions and so drive the structure naturally to this state. 
Morever, once a foam has solidified its structure is set, 
while granular assemblies are fragile by the very nature 
of this state. This difference means that most measure- 
ments in granular assemblies are plagued by large fluctu- 
ations, while foams are amenable to more conclusive ex- 
perimental observations. Therefore, this mapping creates 
a formal gateway between the two systems and makes it 
possible to carry out experiment on one type of systems 
in order to learn about the other. 
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